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Spiral structure is one of the most common structures in the nature flows. A general 
steady spiral solution of incompressible inviscid axisymmetric flow was obtained ana- 
lytically by applying separation of variables to the 3D Euler equations. The solution, 
depending on 3 parameters, describes the spiral path of the fluid material element on 
the Bernoulli surface, whereas some new exact solutions were obtained to be bounded 
within the whole region. The first one is a continued typhoon-like vortex solution, where 
there are two intrinsic length scales. One is the radius of maximum circular velocity 
the other is the radius of the vortex kernel rk = V^rm ■ The second one is a multi-planar 
solution, periodically in z-coordinate. Within each layer, the solution is a umbrella vortex 
similar to the first one. The third one is also a multi-planar solution in z-coordinate. In 
each layer, it is a combination of two independent solutions like the Rankine vortex, which 
is also finite but discontinued for either vertical or horizontal velocity. The fourth one 
is a multi-paraboloid vortex solution finite for z-coordinate but infinite for r-coordinate. 
Besides, some classical simple solutions (Rankine vortex, Batchelor vortex. Hill spherical 
vortex, etc.) were also obtained. The above explicit solutions can be applied to study the 
radial structure of the typhoon, tornados and mesoscale eddies. Both the solutions and 
approaches used here could also be applied to other complex fiows by the Navier-Stokes 
equations. 



1. Introduction 

The exact solutions of the Euler equations for inviscid flow are quite important for 
understanding how the real fluid will flow. But it is a big problem for solving the Euler 
equations due to the nonlinearity. Within 2D (two-dimensional) context, the complex 
potential can effectively solve the irrotational fiow by turning the Euler equations to linear 
ones. For the 2D inviscid flow with vorticity, the vor ticity-stream for mulation can simplify 
the problem to a generalized Helmholtz equation (jBatcheloiill967l). Co nsequently, some 
special explicit solutions were obtained after this simplificat i on (jWang 1991: Saffnia;^ 
19921: iMaida fc Bertozzill2002t IWu et aLll2006l ). And lLou etldl (|2007t) also obtained some 
non-steady solutions from the steady solutions by the general symmetry group theorem. 

Within 3D (three-dimensional) flow context, the problem becomes more complex, es- 
pecially for the 3D flow with vorticity and swirling. The nonlinear advection term make 
the equations hardly to be solved. If the nonlinear advection term vanishes everywhere at 
certain conditions, the fiow is called the Beltrami fiow. Thus the Euler equations could be 
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simplified to linear equations, and some exact solutions might be obtained. The known 
lin iited exact so l utions of the Euler and N avier-Stokes equations were reviewed in detail 
bv lWand (flQOl : IShtern fc HussainI ()l999l) . Such solut ons from simplicity to complexity 
include concentric flows , para llel flows, linear flows, the generalized Beltrami flows in ax- 
isymmetric case (,,Wangj|l99l[ ). Among them, the exact solutions for generalized Beltrami 
flows are the most complex, and they describe swirling flows of widespread technological 
importance. 

Within 3D axisymmetric context, the simplest Beltrami flow is a steady swirling flow, 
where both the radial velocity and azimut hal velocities are de pendent of the radial co- 
ordinate 2, but of the axial coordinate r (iShtern et aZ.lll997l ). In this case, the Euler 



equations for the radial velocity, the azimuthal velocity and swirl are decoupled, where 
the coupling of three velocity components are arbitrarily constructed rather than de- 
rived from reasonable b oundary conditions. Th e solutions include the Burgers vortex, 
the S ullivan vortex (e.g. Wand 199ll : iWu et aLl[2006i ). and a family of general sink vor- 



tices (jShtern et a^.l 119971 ). Another simple axisymmetrically Beltrami flow is symmetric 



sink without swirling. Some of them are impingement of two rotating streams, rotational 



stagn a tion flow over a plate, the Hill spherical vortex and etc ( Wang. .199li IWu et al 



2006[) . iMoffattl ()l969l ) also extended the Hill spherical vortex to allow for a swirl with 



nonzero circumferential velocity that causes a circulation around the z-axis. 

Although the above solutions are helpful, most of the symmetric sink solutions are 
decoupled and unbounded, except for, e.g., the Hill spherical vortex and its extension. 
The decoupled and unbounded solutions mi ght be of some limited physical interest. To 
avoid this arbitrariness and unboundedness. Long (jT958, 1961) tried first to obtain the 
coupled equations for conical similarity swirling vortex, and the Euler equ ations can 
be simplified to the Bragg-Hawthornc equation, or the Long-Squire equation ([Batchelor 
[1967; Saffman 1992; Frewer et a l. 2007). However, such equation is generally nonlinear, 
no exact explicit soluti on has been obtaine d. 

In contrast to that, iFrewer et al] (|2007l ) used Lie-algebras to investigate the Bragg- 
Hawthorne equation. They found out eight distinct local Bragg-Hawthorne symmetries 
within five algebras. However, such symmetries ar e different and even more symmetries 



than the Euler equations do (jFrewer et al\\2001\ ). To dispel the inequivalent between 



them, they need to obtain some additional equations. Even though, th ey also find only 
the simple linear Bragg-Hawthorne equation can be solved explicitly (|Batchelorlll967 : 
Frewer et allhOOlh . 



As an alternative way to avoid the above disadvantage, we conside r ed to solve the cou- 
pled nonlinear Euler equations following the deviation by iBatcheloi ( 1967 ) , which is the 
most complex Beltrami flow with that all the velocity components are dependent on both 
the radial coordinate r and the axial coordinate z. Due to the com plexity, the general 
expli c it solution for s uch flow still lacks even in the recent literatures (jMaida fc Bertozzi 
2OO2L Wu et al. 2006 ). In th is study, we solved the Euler equations by applying separa- 



tion of variables ( Sun 201 1[ ). As a result, a general exact spiral solution was presented 
in §[2l some special explicit solutions were given in §[3l Discussion and conclusion were 
respectively given in §|4] and in §[5l 



2. General solution 

We considered the steady solution of the incompressible Euler equations for axisym- 
metric flow at present study. It is convenient to use a cylindrical coordinate system 
{r,9,z) with the velocity components (K-jV^, Vj), and all the velocity components are 
the functions of r and z but 9, due to the axisymmetric. As Vr = Vr(r, z), Vg = Ve(r, z) 
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and Vz — Vz{r, z), the governing equations, including mass-conservation and momentum 
equations, are 

d{rVr) d{rVz) 



= (2.1a) 



dr dz 
^ dr ^ dz r 



or oz r p or 

Or Oz p oz 

It should be noted that there is no length scale in Eq. (|2.1l) . t hus the so l ution can be 
uniformly stretched by simply multiplying a complex constant. iBatcheloii (|1967[ ) noted 
that Eq. p.lb .) and Eq. ()2.1b ) can be solved through the introduction of the axisymmetric 
stream function '^{r,z). We tried to find the solution of the above equations by separation 
of variables \I'(r, z) ~ R{r)H{z). One such solution can be written as, 

K = ^H\zl Ve ^ X^Hizl Vz ^ -^Hiz) (2.2) 
r r r 

where ' presents first deviation and A is a complex constant. In this way, equation (j2.1b .) 
and equation (|2.1b ) are satisfied automatically. The path of a fiuid material element can 
be obtained, 

lnr(0) = (2.3a) 

^lir^z) ^ const. (2.36) 

In r — 6* plan, it is a logarithmic spiral {H' ^ 0), except for ' = (a circle). So we called 
the solution is spiral solution. In fact, the path is right on a Bernoulli surface given by 
Eq. (|2.3b ). on which the streamfunction is defined (Batcheloriil967i) . 
Equation p. lb ) and equation (|2.1d ) become, 

(^)(:^)'i7'2 _ i\^)HH" - %R^H^ = -^^ (2.4a) 



r 



j,3 p Qj. 



i^fHH'-i^)i^yHH' = -'-^ (2.46) 
r r r p oz 



Hence, the pressure can be solved from Eq. (|2.4b ). 

p(r, z) = lpH'[i^)i^y - i^f] - pQ{r) (2.5) 
2 r r r 

Substituting Eq. (|2?5l) into Eq.(|2^). it yields 

(^)(^)'if'2 _ (^)(:^)(^^.) _ i^R^)H' hi^r - {\^y]'H' + Q'{r) (2.6) 
r r r r r' L r r r 

Recall that R and H are independent, we can solve the above equation. There are three 
kind of functions for li[z) in Ea. (l2.6p . The first trivial one is iJ = 1, any differentiable 
function R(r) would be the solution, which we discussed in S I3.2.1I Besides, it also yields 
two kind of non-trivial solutions for H{z): a) both H'^ and HH" are independent to H^, 
and b) both H'^ and HH" are proportion to H^. This yields 

HH" = —^H'^,n = 1, 2, or, (2.7a) 
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HH" = k^H^,and, H'^ = k'^H^ ± AajSk^ . (2.76) 
where k and P are two complex constants. The solutions are obtained for H{z) and Q{r), 

H{z) - z", Q{r) = h^)^ ~{n~l) [ ^dr, or, (2.8a) 
2 r J f-^ 

H{z) = ae'^^ ± ;3e-'=^ g(r) = ±2aPk^{-f, k^ ^ 0, (2.86) 

r 

Substitution Eq. (|2.8p into Eq. (|2.6p to eliminate the function for z, and letting 7 = A^+fc^, 
it yields 

-5^' = [(v)'-(7)(v)T (2-9) 

where yfc = is for Eq.dl^) and fc 7^ is for Eq.dHJj). The solution of Eq. dH]) for 
7^0 yields, 

E(r) = ar^e"^'^'"' (2.10) 

and the solution for 7 = yields 

R{r) ^ar^ ~b,or (2.11a) 
R{r)^ae~^''' (2.116) 

where a, 6 and c are complex constants. And the solution Eq. (j2.11b ) is independent to 
Eq. (|2.11b ') under the condition of 7 = 0, as the function Eq. (|2.9p is nonlinear. Accord- 
ingly, by integrating Ea. (|2.5p with Eq. (|2.10p the pressure is, 

p{r,z)^ pH^{d^4R^)^ pQ{r) (2.12) 

Thus, Eq. f^ with Eq. ((^ and Ea. (pTTU|) or Eq. ipH]) give the exact solutions of the 
flow velocity. It is obvious that there are three independent parameters (e.g., a, fc. A) in 
the solutions. And other new solutions can also be obtained by combining the different 
solution at different regions, like that of the Rankine vortex. 



3. Special solutions 

As we know, the Euler equations are defined in the real field. While the above general 
solutions are defined in the complex field. If the parameters in the solutions are real, 
the solutions are real. If the parameters are complex, we should try to take the real 
part of the complex as the solution. However, such approach can not work in general 
(some additional conditions are required), in that the Euler equations are nonlinear. 
Such kind of property is common for nonlinear equations, even for a sim ple nonlinear 
algebr a equation. This problem is similar to the previously mentioned one in lFrewer et ali 
(|2007l ). who also found out not all the Lie-algebra solutions satisfy the Euler equations. 
In their studies, some additional equations required. 

For the velocity, the combination of Eq. p.8p and Eq. p.lOp or Eq. (|2.1ip according to 
Eq. (|2.2p gives the solution. However, there are too many parameters (including A, a, 
/3, 7, fc, a, 6, c, n) in such solution. Our further investigations also pointed out that 
Eq. l|2.1b ) holds only for A being a real number, i.e., Xi — 0. Hence, we restricted our 
following investigations toa = l,/3 = 0,a = l and 6 = without loss of generality. 
Although such restrictions make the flows to be much more simple, they are not trivial. 
In the following section, we gave some interesting solutions to show how the fluid flows 
in such conditions. 
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Figure 1 . Paths of the fluid material elements for the umbrella vortex. 

3.1. k^Q 

For fc 7^ 0, H{z) = e^^ , ■^r = — kf + - Af and 7^ = 2{krki + A^A^), we have the 
parameters 

1 



— {arki + ttikr), Bi = (a^-Xi + a^Ar), Ci — Qi - 
And the solution is 

Vr — {A,. COS (f) ~ Ai sin (j))re'^ , 
Vg — [Br cos (f) — Bi&u\(j))re^ , 
= -2(C,. cos - C, sin 0)e'^, 



Or - -{ar-fr - aiji)r 



1 



(3.1a) 
(3.16) 



(3.2a) 
(3.26) 
(3.2c) 



where = kr 



and 6 — fc, 



^7,;r^. Although either fc or 7 might be complex 



number, only some especial combinations could be the solutions. The first two types 
require 7i = 0, in one case the magnitude of the circular velocity Ve can be either less 
or larger that of the radial velocity K, and in another case the magnitude of Vg is larger 
than that of Vr- Both suit for typhoon structure. The last two types require 7^ = 0, 
where the magnitude of Vg is less than that of Vr- 

3.1.1. kr 7^ 0, ki = and 7^ = 

In this case 7 is a real number with 7^ > 0, and fc^ = implies a exponential changes 



V7. — k^rG 
Vg = Xrre'' 

V, = -2(1 



1 



;lrr )e 



(3.3a) 
(3.36) 
(3.3c) 



The solution is finite within the whole domain, if we choose appropriate fc^ > for z < 
and fcr < for z > 0. The circular velocity Vg is the same with that of the Taylor vortex 
for any fixed z. The circular velocity maximum is Vg„i — XrVmc''^^^^^'^ at = 4/7^. 
And the vertical velocity Vz changes it direction at = 8/7^, which we call as the 
radius of vortex ker nel . Such vertical velocity dist ribution is something like that of 
the Sullivan vortex ( Tone et"allll994l : Wu et al. 20061 ) . The paths of the fluid material 
elements are illuminated in Fig [TJ In r — plane, the paths are cyclonic logarithmic 
spirals (Fig[T^) or anticyclonic logarithmic spirals (FiglTja) with Inr = krO/Xr . In 2: — r 
plan (Fig[T};), the fluid material element moves spirally along the surface decided by 
z = {■^'Jr'r'^ — 2\nr)/kr + const. Figure [Hi also shows the path in the 3D space. In the 
inner region r < rk, the flow cyclonically ascends from below to above as shown by 




upward arrow. In the outer region r > rk, the descends apart from the kernel. So the 
mass conserves in this case. 

^e'^'-^ = - / 2'KrV^dr = / 2T:rV^dr (3.4) 

e Jo Jrk 

There are two intrinsic length scale in the solution. One is the radius of maximum 
circular velocity r„j, the other is the radius of the vortex kernel — \/2rm- According 
to the shape of the paths, we called it the umbrella vortex. This new vortex solution is 
very like that of the typhoon structure, which we know quite limited about. 

3.1.2. kr — Q, ki ^ Q and 7,; — 

In this case 7 is a real number with 7^ > (A^ > kf), and fc^ = implies a sinusoid 
function for z by taking fc^ > without loss of generality, 

Vr ^ -kism{kiz)re~^^'^''^ , (3.5a) 
Ve = XrCos{k,z)re-^'^'-'^\ (3.56) 

V, = -2(1 - l-/rr^) cos{hz)e-i"''-''\ (3.5c) 



Similar to the above solution, Eq. p.5p is finite within the whole domain. But the present 
solution has multiply layers by noting that the solution is periodic in z-coordinate. The 
fluid material elements are restricted within different vertical layers. So we call such flow 
as multi-planar flow. In each layer (e.g., ^ kiZ ^ ^), the flow has a similar behavior 
like that in Ea. (l3.3p . except for that there are both inflow {hz > 0) and outflow {kiZ < 0) 
at present solution. A similar flow path of the fluid element can be find in Fig[T] There 
are also two intrinsic length scale in the solution. One is the radius of maximum circular 
velocity r,n, the other is the radius of the vortex kernel rk — V2r,n- Another notable 
property of the multi-planar solution is that there might be numberless solutions for a 
given boundary condition at z-coordinate. So the solution of the 3D Euler equations for 
a given boundary condition might not be uniqueness. 

3.1.3. kr — 0, ki — ±Xr and 7^ = 

In this case, 7^ = 7i = 0, the magnitude of Vr is the same with that of Vg (like Ekman 
spiral). And the solution of R is Eq. (l2.11l) for 7 = 0. In Eg. (12. lib ), the solution is finite 
for inner kernel r < but infinite for far region r — >■ 00. While in Eq. (j2.11b ). the solution 
is finite for far region r > but infinite for r — > 0. Similar to the solution of the Rankine 
vortex, we can combine the above solutions to obtain a uniform solution at the interface 
r = Tj, which is finite at whole region. However, either radial velocity Vr or the vertical 
velocity Vz might be discontinued in Eq. (|3.6|) . Additional notation for rk will be given in 

§[3221 

For example, we choose K to be continued at the interface rf = e^'^^' = 1/c, hence 
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Figure 3. Paths of the fluid material elements for multi-paraboloid vortex. 



the solution is, 



Vr — kir sm{kiz), Ve — ±kir cos{kiz), Vz — 2cos{kiz) 



(3.6a) 



Vr = —e-"''' sin(fciz), Vg ^ ±— e""'"' cos(fc,z), V, = -2ce-"''' cos(fc,z). (3.66) 
r r 



The present solution has same function as Eq. (|3.5p in z-coordinate. Besides, the velocities 
change their directions along the r-coordinate, which is like that in Ea. p.51) . Figured] 
shows the path of a fluid material element within — 7r/2 ^ kiZ ^ tt/2 for Eq. (|3.6b ). 
According to Eq. (l2.3p . in r — plan, the paths are logarithmic spirals Inr = ta,ii{kiz)0 
(Fig[5^,b). In z — r plan (Figl^t), the fluid material element moves spirally along the 
surface decided by cos(fciz)r'^ = const. Figure also shows the path in the 3D space. 
It should be noted that the flow field is respectively inflow and outflow at lower part 
— 7r/2 ^ kiZ ^ (Fig[2^) and higher part ^ kiZ ^ 7r/2 (Figl^b)- However, the path of 
the fluid material element is always cyclonic from below to above (Figl2}l). 

The radial velocity Vr and the circular velocity Ve continue, but the vertical velocity 
Vz discontinue at r = r^. If the flow within the inner kernel ascends at r < r^, then the 
flow descends at r > r^. The total updraft mass is equal to the total downdraft mass. 



2'Krf cos{kiz) 



2'KrVzdr 



2nrVzdr 



(3.7) 



3.1.4. kr ^0, ki^O and 7,. = 

In this case 7 is a pure image number, and 7^ = implies kf. — kf + = 0. The 
solution is. 



Vr = [kr cos{hz - ^Jir'^) - h sm{kiZ - ^■^ir'^)]re'''-\ 
Ve = A,. cos{hz - \^ir^)re^-'-'' , 

o 

Vz = -2[cos{kiZ ~ --lir'^) + —r^ sin(fc,z - -"f,r^)]e^-\ 



(3.8a) 
(3.86) 
(3.8c) 



It is obvious that this solution can only be an inner solution r < rk, as they might be 
infinity for r — )• cx). In this case, the magnitude of Vr is ^/k^+kf, which is larger than 
that of Ve- So the rotation is relatively small, which can be seen from Fig|3^,b. Both 
the radio velocity and azimuth velocity may change directions along the r-coordination 
and z-coordination. The solution is also a multi-paraboloid vortex by noting that the 
velocities are spatially periodic shift along the paraboloid of kiZ = \'yir'^ (Figllt). For 
example, Fig|3ji depicts the the paths of the fluid material elements within a paraboloid 
layer. Due to the smaller azimuth velocity, the paths only have some slight twists, whereas 
there are lots of spirals in FigHJ 
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3.2. k = Q 

In this case, let H{z) — z",n — 0, 1,2. As the solution may be infinity for n = 1,2 as 
z CO, so it should be taken only as a solution in lower layer ^ < oo. 

3.2.1. H{z) ^ 1 

Firstly, considering n = and H(z) — 1, the radial velocity Vr vanishes according to 
Eq. (|2.2l) . This is barotropic flows in the geophysics, as the solution is free of z-coordinate. 
Any differentiable function R{r)^ the circular velocity Ve and vertical velocity Vz would 
be the solution from Eq. (|2.6p due to that Ve and Vz are decoupled for H{z) = 1, which is 



known as a stretch-free inviscid vortex can have arbitrary radial dependence (|Wu et al 
[2OO6.I . In present solution, Ve and Vz does not fully decoupled, as Eq. p.Qj) shows. 

R{r) _ R'{r) 



Vr - 0, = A- 



Vz - 



(3.9) 



Given R{r) = 1 is the line vortex, R{r) = is the solid rotation. And the Rankine vortex 
and the Taylor vortex can be obtained by given R{r) = ar^ + b and R{r) = r^e"'"' , 
respectively. The Bat chelor vortex an d Oseen-Lamb vortex can also be obtained by given 
R{r) = 1 - etc (|Wu et oi.ll20(ih . 

Although any differentiable function R{r) gives the vortex solution, not all the vortex 
can be hol d due to the instabilit i es, including th e centrifuga l in stabil i ty an d the shear 
instability (jCriminale et a/.ll2003t IWu et a/.|[2006t [Su3l2006l) . In[Su3 (|2006l ). a general 
stability criterion was obtained for such circular flows. It was found that the longwaves 
(e.g. wavenumbers are 1, 2, etc.) might be unstable and break into asymmetrical vortices 
if there is a local maximum in vorticity distribution along the radial coordinate. Th e 
numerical simulations also approved this ( Roger et al. 199Ct Belotserkovskii et al. 20091 ). 



3.2.2. H{z) ^ z 

Secondly, consider 71 = 1 and H{z) = z. If A 7^ 0, the solution is, 

1 



Vr 



are 



Ve = aXre 



z, Vz = -2a(l- -AV^)e 



(3.10) 



And if A = 0, the circular velocity Ve vanish due to this. Similar to Eq. p.6|) . the present 
solution Eg. p. Ill) is also the combination of two parts at r = r^, where a discontinuity 
might occur. One should note that we have chosen a — b — I and c — 1/r'j in Eg. p.Gp . 
Hence, the solution for for the whole domain is. 



Vr 



Vr = ar, Ve =0, Vz = -2az 
-- -e""''', Ve =0, Vz^lbcze: 



(3.11a) 
(3.11&) 



The vertical velocity Vz discontinue at r = r^, and the vertical velocity discontinuity is 
AVz = 2a(l -I- brf). Similarly, we can choose a — — fece^"^''*, thus the vertical velocity is 
continued but the radial velocity has a discontinuity of /Wr = a{c/ri + Ti), which is also 
a two-cell vortex. 



3.2.3. H{z) = z2 

Finally, consider n — 2 and H{z) 
Ve vanish due to this. Only R 



^. As 7 = implies A = 0, the circular velocity 
could satisfy Ea. (l2.6p . and the solution is, 

Vr = 2rz, Ve = 0, 14 = -2^^ (3.12) 

This is the soluti on of rota tional s tagnation flow over a p late (' Wandll99lh . or the Hill 



spherical vortex (|Batchelor..l967; .Tone et a/.lll994l: IWu et qL, 20061) . 
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4. Discussion 

4.1. Extensions of Solutions 

In above section, we assumed & = in Ea. (|2.11b ) to obtain the inner solution. Alterna- 
tively, we can take a = and 6 ^ to as a outer solu tion, which is a well-known pure 
vortex solution for stretch-free vortex dWu" et aL|[2006l ). It is obvious that Eq. (j2.11b ) is 
the solution of the Couette- Taylor flow, we can also use Ea. (|2.11b ,) and H{z) to obtain 
new solutions. 

In Eq. (|3.5|) and Eq. (|3.6p . both solutions have a same ^-coordinate dependence function, 
so these solutions can also be combined for some new solutions. For example, if we take 
Eq. (|3.5l) for r < ri and Eg. p. 6b ) for r > as a new solution, where rf = 1/(^ — 1) 
and c — {1 ~ + ^) in Eg. (13.6b '). both the Vr and dVr/dr of the new combined 

solution are continued at r = r^. As this new solution is continued, it might be better 
than that of discontinued solution in Eq. (|3.6|) . It is noted that the value of 7r could be 
negative in the near- kernel region for present solution, although the solution of Eq. ()3.5p 
requires 7^ > in the far region. The negative 7^ implies a very fast tangential velocity 
rise in the near-kernel region, which is true in the typhoon. 

According to the previous studies ( Batchelorlll967tlFrewer et~al. 20071) . the solu t ion se t 



is quite large, which can also be seen from the Ea. (|2.2p . According to iBatcheloij (|l967l ) 



the vertical velocity Vg = C(4')/r, while we simply took the velocity components as a 
form Ve = AvE'/r in Eq. (j2.2p . Beyond present study, there should be other axisymmetric 
solutions. For example, we can apply the same approach to find other exact solutions 
by taking Vg = X/r or Vg — A4'^(r, z)/r, etc. The method used in present work could 
also be applied to other complex flows, e.g., the geophysical fiow in a rotating frame 
(/-plane), the mag netohydrodyn amics (MHD) in astrophysics, and even for the viscous 



flows. According to IWangj (|1991l ). the generalized Beltrami flow might also be the steady 
solution of Navier-Stokes equations, if the curl of vorticity is potential. Moreover, the 
present solutions might further be used to obtain non-steady solutions of Navier-Stokes 
equations, like that of Oseen-Lamb vortex. 

4.2. Possible Applications 
As mentioned above, a stretch-free inviscid vortex (two-dimens ional axisymme tric colum- 



As mentioned above, a stretcn-tree mviscicl vortex (two-dimens ional axisymme tric colum- 
nar vortex) can have arbitrary radial dependence iov H{z) — 1 ( Wu et allbOOGl ). However 



many well-known vortex solutions are always similar to either Eg. (12. lib ) or Eq. (|2.11b ). 
It is from this study that these two-dimensional axisymmetric columnar vortices are also 
the three-dimensional axisymmetric columnar vortex solutions. So we can find either the 
Rankine or Taylor vortex for any fixed layer z = const. 

As some new exact solutions are finite within the whole region, they can be applied 
to study the radial structure of the typhoons, the tornados and the mesoscale eddies 
in the geophysical flows. As mentioned above, there are 3 independent parameters in 
the solution. To determinate them, we need the measurable qualities, such as the radius 
of circular velocity maximum rm, the maximum circular velocity Vsmj the spiral angle 
defined by k/X, the descend/ ascend flow mass in Eq. p.4|) and Eq. p.7p . the horizontal 
infiow/outfiow mass (which can also be determined by both k/X and the descend/ ascend 
flow mass), the angular momentum m, etc. According to Eq. (|2.2p . the angular momentum 
m = rVe = -XR{r)H{z), thus the solutions of R{r) (Eq. ipiB and Eg. dCTl ) can be 
applied to discuss the angular momentum of the vortex. For the typhoon observations, 
such solutions can be used to fit the real velocity distribution along the radial coordinate. 
This may also be useful to classify the typhoons and the tornados according the flow 
structures provided by above solutions. 
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5. Conclusion 

A general exact axisymmetric spiral solution was obtained for 3D incompressible Euler 
equations, and some exact solutions were also obtained. The solution describes the spiral 
path of the fluid material element on the Bernoulli surface. The explicit umbrella vortex 
and multi-layer vortex solutions, which are new in this investigation, might be used to 
describe the 3D structure of the tropical cyclones, tornados and mesoscale eddies in 
the geophysical flows. The solutions also imply that the spiral structure is the intrinsic 
structure of the flows in the nature. 

The method used in present work is very straightforward, and it could also be applied 
to other complex flows, e.g., the geophysical flows in a rotating frame (/-plane), and even 
for the non-steady viscous flows. 
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problems with the author and encouraged the author to finish this work. The author 
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acknowledged, who led the author to this field. This work is supported by the National 
Basic Research Program of China (No. 2Q07CB816004), and the Knowledge Innovation 
Program of the Chinese Academy of Sciences (Nos. KZCX2-YW-QN514). 
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